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Gamma-ray bursts (GRBs) have been invoked to explain both the 511 keV emission from the 
galactic bulge and the high-energy positron excess inferred from the ATIC, PAMELA, and Fermi 
data. While independent explanations can be responsible for these phenomena, we explore the 
possibility of their common GRB-related origin by modeling the GRB distribution and estimating 
the rates. For an expected Milky Way long GRB rate, neither of the two signals is generic; the 
local excess requires a 2% coincidence while the signal from the galactic center requires a 20% 
coincidence with respect to the timing of the latest GRB. The simultaneous explanation requires a 
0.4% coincidence. Considering the large number of statistical "trials" created by multiple searches 
for new physics, the coincidences of a few per cent cannot be dismissed as unlikely. Alternatively, 
both phenomena can be explained by GRB if the galactic rate is higher than expected. We also 
show that a similar result is difficult to obtain assuming a simplified short GRB distribution. 



Gamma ray bursts are among the most energetic 
events in the universe. Although there is a significant 
correlation of long GRBs with star-forming metal-poor 
galaxies [l], many long GRBs are observed in high- 
metallicity galaxies as well Therefore, one expects 

that long GRBs should occur in the Milky Way Galaxy 
once every 10 4 — 10 6 years [rjl-flol|. This large uncertainty 
in the rate is largely due to selection effect, as the ob- 
served GRB rate and the opening solid angle necessary 
to determine the absolute rate are sampled from two dif- 
ferent populations. Furthermore, the large variations ob- 
served in the beaming angle, and the uncertainties in the 
measurement of those angles make it difficult to accu- 
rately constraint the absolute GRB rate in the Milky Way 
7, 10]. Short GRBs, which probably result from com- 
pact star mergers, should occur in every galaxy, includ- 
ing the Milky Way, at a rate that is comparable, although 
slightly lower. Short GRBs are likely less beamed then 
long GRBs, and some even appear to not be beamed at 
all [Ill4l3j . For the purpose of this work, we follow limits 
set using observation of > 0.1 mJy sources at 8.44 GHz by 
the Very Large Array (VLA) , and we assume a lower limit 
on the beaming angle of a few degrees [3] . This sets an 
approximate upper limit on the rate i?GRB ^$ 0.5 x 10~ 4 
per year. 

GRBs shoul d p roduce high energy electrons and 
positrons [8l ll5lll6j . and it is therefore natural to consider 
GRBs as a possible explanation of recent astro-particle 
observations involving electrons and positrons [1, [Tt], EH , 
as well as cosmic rays (l9j . 

High energy electrons and positrons are produced dur- 
ing the initial burst [2(J[2lJ]. However, due to the short re- 
combination time scale, these electrons are unlikely to es- 
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cape the initial fireball [l8j]. On the other hand, high en- 
ergy electrons and positrons can also be produced by the 
interaction of the 7-rays with the surrounding medium 
in the jets [H, [Hj], or ahead of the fireball fia liH |24| . 
In the former case the annihilation of the pairs freezes 
out as the jet expands, and electrons and positrons are 
allowed to escape to the outside. In the latter case 7-ray 
photons are back scattered into the ejecta and allowed to 
pair produce through interactions with subsequent GRB 
photons. In both cases conversion efficiencies is expected 
to be on the order of 1 — 10%. 

Results from PAMELA have shown an excess in the 
positron fraction above 10 GeV [251 ] ■ while Fermi, PPB- 
BETS and ATIC have shown an excess in the total elec- 
tron and positron flux above 100 GeV 26 28]. Both of 
these results can simultaneously be explained by a nearby 
(~ lkpc) GRB-like event 2 x 10 5 years ago (!§]. Most 
GRB models usually predict the electron positron pairs 
to form with MeV energies [HI, [24J. However, the in- 
teraction of backscattered TeV photons created during 
the burst with optical photons from the afterglow would 
also give rise to a GeV — TeV population of electron and 
positrons. The total energy of the produced pairs would 
be comparable to that of the initial high energy photons: 
10 50 erg for a typical spectral index of 7 = 2.2 [291 ]. 

Furthermore, observations by INTEGRAL satellite of 
the central region of the Galaxy have discovered a spec- 
tral line of energy 511 keV coming from the galactic 
bulge [30M35| . This emission could also be explained by 
GRBs inside or near the Central Molecular Zone (CMZ) 
occurring every ~ 10 4 years [UHl]. Indeed, a conver- 
sion efficiency on the order of 1 — 5%, would generate a 
number of MeV positrons 




where E is the total initial energy of the burst fl5Hl~8l l2~i| . 
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The electrons and positrons lose energy via ionization on 
a time scale of ~ 10 7 years and diffuse through the bulge, 
traveling a distance similar to the extent of the 511 keV 
image before stopping [T3)I3- Thus, the explanation is 
in agreement with the morphology of the observed signal. 
Since the probability of annihilation is much higher at low 
energies, most of the positrons injected into the central 
region slow down and annihilate at rest [171 Il8| ] . 

The explanations put forth to explain these two phe- 
nomena in Refs. fPH, Il8| and in Ref. [2!| have not made 
any contradictory assumptions. The morphology and the 
intensity of the 511 keV line depends on some (reason- 
able) assumptions regarding the gas density and mag- 
netic fields in the Galactic Center [l7l . fl8l ] , which do not 
affect the high-energy signals [2{|. However, it has not 
been shown that the rates of GRBs and their distribution 
in the Galaxy afford simultaneous explanation of both 
the low-energy and the high-energy anomalies. This is 
the question that we will address. 

Gamma ray bursts form to distinct populations, long 
and short 36]. Although there remain a great deal of 
uncertainties regarding the origin of GRBs, long GRBs 
are usually associated with core-collapse supernovae in 
star forming galaxies [37J, while short GRBs are usually 
attributed to the merger of neutron stars and/or black 
holes [H| . We will explore the possibility that both type 
of GRBs may be responsible for the 511 keV and the ex- 
cess measured by ATIC, Pamela and Fermi, simultane- 
ously. 

We model the history of Galactic GRBs in a Monte 
Carlo simulation by randomly generating bursts to de- 
termine the probability that both the electron-positron 
excess in the ISM, and the 511 keV line might originate 
from GRBs. To this end, we calculate the probability 
that a measurement made at a random time in the his- 
tory of our galaxy falls within the required time con- 
straints. We will assume that the spatial distribution of 
long GRBs follows the distribution of stars determined 
from the star counts [3!l|4(|. For the short GRBs, which 
occur farther in the halo, we will use the distribution of 
sources obtained from observations [4l|. We show the 
simulated maps of the GRBs locations in FIG. [TJ 

It is reasonable to assume that the distribution of long 
GRBs in the galaxy is closely related to the stellar density 
profile. We model the stellar distribution in the galaxy as 
a two-components system comprised of a thin disk and 
a bulge. Below we will also include the effects of the 
spiral arms. We neglect the spheroid component because 
it accounts for only 1 — 5% of the mass of the galaxy. 
We use the stellar density profiles of the Bahcall-Soneira 
model HEl: 



n D = Pn(ro) / drdzdip I dM $(M) x 



re - z /H(M) e -(r-r )/h 



(1) 

n B = p B J dRdOdip J dM<f>{M)R 02 sin 9e' R \ (2) 
Where z is the perpendicular distance from the plane of 
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FIG. 1: Simulated map of GRB locations in galactic coordi- 
nates for a distribution that follows the star counts [39l . l4p| ] 
(top), and for a distribution inferred from observations of 
short GRBs HH (bottom). 



the disk, r is the galactocentric distance with tq = 8 kpc 
the galactocentric distance of the sun, and R is the spher- 
ical distance from the center of the gala xy. Following the 
literature we take h — 3.5 kpc [39l l40j|. The normaliza- 
tion constants po and pg are only relevant in our work 
for determining the relative fraction of stars in each com- 
ponents of the Milky Way. Based on mass studies, we 
assume that 85% of the stars are located in the disk and 
the remainder in the bulge. The luminosity function is 
given by 



$(M) = 



n,10 g ( M - M ') 

[i+io 5 ( q -' 3 '( m -- m *)] 1/ ' 5 ' 

const = $(15), 

0, 



-6 < M < 15 

15 < M < 19 

M < -6 or M > 19 



where n, = 4.03 x 1(T 3 , M* = 1.28, a = 0.74, f3 = 0.04, 
and 6 = 1/3.4. Finally, to facilitate the random sam- 
pling process, we adopt a simplified version of the scale 
height. Old stars with visual magnitude M > 5.1 are 
generally associated with an exponential scale height 
of 325 pc. Main sequence stars with visual magnitude 
M < 2.3, are well described by an exponential scale 
height of H = 90 pc; for 2.3 < M < 5.1, the scale height 
is interpolated with a linear regressions [1^, H^j . Instead 
of performing this regression we will, for the purpose of 
this work, assume the scale height to be well described 
by a step function. In other words: 
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H(M) 



90 pc M < 3.7 
325 pc M > 3.7. 



We will model the distribution of short GRBs based on 
a simplified star formation disk model [42] . This assump- 
tion, although not completely accurate, reflects the fact 
that a significant fraction of short GRBs could originate 
from massive star collapse [43J, |44| . 



n Sh = Psh / dRdOdip Re 



-1.670(H) 



(3) 



where 0(R) = j^— , with i?haif the half-light radius of the 
Milky Way. We treat each coordinate as completely inde- 
pendent and, after normalizing the stellar density profiles 
such that the area under the curves is equal to one, we 
treat the distributions as probability density functions. 
For the purpose of the Monte Carlo we distribute the 
bursts over time using a Gaussian profile with average 
rate 



R. 



1 



10 4 yr 



(0.1-0.5) 



1 



10 4 yr 



(4) 



We first proceed with the long GRB analysis. We ran- 
domly generated 10 6 bursts based on the spatial and tem- 
poral distribution described above, and we recorded the 
number of events (i) occurring within 1 kpc of the sun, 
as well as those (ii) inside CMZ or beamed towards the 
CMZ. We find that about 0.2% of all events fall in cat- 
egory (i). Similarly, we find that about 12% of all even 
fall in category (ii), and hence, the average rate of GRB, 
and, therefore, of electrons and positrons injection in the 
central part of the bulge is 



R i 



1 



8.5 x 10 4 yr 



(5) 



Let us now discuss the probability of a GRB within 
1 kpc of the sun. For a uniform distribution of stars in 
the disk, we obtain the rate of nearby GRB 



R 



1 



0,uniforn 



5 x 10 6 yr V 1.0 



(6) 



However, the sun is located in a spiral arm, and the 
local density of stars is expected to be somewhat higher. 
Following Ref. [45[ , we model the Milky Way spiral arms 
of the form: 

r(9) =q e (e - e ° )tanl , 

where qo = 2.3 kpc, i = 13.8°, and 6q ranging from to 
in steps of | representing the starting angle of each 
arm. The brightness of the arms varies with the passband 
it is observed in. In the O passband the contribution 
of the arms to the total brightness increases from 17% 
at 3 kpc to 50% at 15 kpc; those contributions are 20% 



stronger in the B passband, and 50% stronger in the A 
passband (4(| . Given the location of the solar system, we 
find that this over-density results in a factor / sp i ra i ~ 2 
enhancement in the local GRB rate: 



R, 



/spiral \ 



R 



1 



0, uniform 



2.5 x 10 6 yr 



( r 4,GRI 

V 1.0 



(7) 

Let us now compare the rates in Eqs. ([S]) and J7]) with 
the mode requirements. As discussed in Ref. [18j, the 
positron population in the bulge enters a steady-state 
regime at the right density to explain the 511 keV for 
R CMZ > (3 x 10 4 yr)~ 4 . For a lower rate, the positron 
population rises and dies off after each GRB. Matching 
this rate to Eq. ([5]) requires r 4 GRB « 3, which outside the 
expected range r 4 GRB =0.1 — 0.5 Q. 

Of course, there is no reason to believe that the rate of 
511 keV photons is time-independent and that we don't 
just happen to observe it at the peak shortly after a GRB. 
In this case, one can estimate the probability of such a 
coincidence. If GRBs are responsible for the 511 keV line, 
and if one takes r 4 GRB =0.5, than one requires that the 
current observation be within the last 20% of the average 
GRB time separation, which is about a 20% coincidence. 

For a local GRB to explain the high-energy electrons 
and positrons, a GRB should have happened in the solar 
neighborhood about ~2x 10 5 yr ago [29(. To match the 
rate in Eq. 0, one must assume r 4 GRB sa 12, which is, 
again, outside the expected range r 4GRB = 0.1 — 0.5 0]. 
We estimate that a ~ 2% coincidence is required for a 
sufficiently recent GRB to explain ATIC PAMELA and 
Fermi data. 

Finally, we ask what would be required to explain both 
observations by GRB. If a GRB occur at a much higher 
rate, i? GRB = 12/10 4 yr, then both observations have a 
natural explanation by GRBs. If the GRB rate is lower, 
as expected, namely i? GRB = (0.1 — 0.5)/10 4 yr, then a 
2% coincidence is required to observe a local GRB suf- 
ficiently recently while the same requirement applied to 
a cenral Galactic GRB requires a 20% coincidence. Of 
course, a simultaneous explanation requires a ~ 0.4% co- 
incidence in this case. This probability is by no means 
small. Many independent observations are carried out in 
search of new physics, including dark matter. The large 
number of independent searches increases the "multiple 
trials" effect, so that one should expect 0.4% effects to be 
promoted by multiple trials. Finally, it is possible that 
the rate mismatch is ameliorate d by the model uncer- 
tainties discussed in Refs. [171 Il8l l29j. 

We repeated the previous procedure for short GRBs, 
ignoring the overdensity effects of the arms, and found: 



j^sh __ _ 
CMZ _ 4.5 x 10 4 yr 



1 / r 4,GRB \ 

10 4 yr V 1.0 / 



R 



■sh 



i 



1.5 x 10 8 yr V 1.0 



(8) 
(9) 



Therefore, in order to explain the 511 keV we need 
r 4 GRB sa 1.5 or a 35% coincidence assuming r 4 GRB = 0.5, 
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while in order to explain the electron and positron ex- 
cess we need r 4 GRB 750, or a 0.067% coincidence. In 
other words, the simultaneous explanation would require 
in this case a 0.023% coincidence. This results implies 
that the short GRBs are less likely to account for both 
effects. 

We conclude that a simultaneous explanation of the 
511 keV line on one hand, and ATIC, PAMELA and 
Fermi on the other hand, based on GRBs in the Milky 
Way Galaxy is possible in one of two ways. First, it is 
possible if the long GRB rate in our Galaxy is 3/10 4 yr, 
which is higher than expected. Second, it is possible if we 



are observing the effects of two very recent GRB, one in 
the Galactic Center and one in the local neighborhood, 
each of which happened more recently than average. The 
probability of the latter is about ~ 2%. It is however un- 
likely that short bursts be simultaneously responsible for 
both effects. Although short GRBs could easily explain 
the 511 keV line, it is hard to reconcile them being the 
cause of the local electron and positron excess. 
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91ER40662 and NASA ATP Grant NNX08AL48G. A.K. 
thanks Aspen Center for Physics for hospitality. 
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